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1. Introduction 

Wigner-Dyson universality of the local correlation of energy levels among various stochastic 
[1] and quantum-chaotic systems [2] under well-defined conditions was established through the 
ten-fold classification of symmetric spaces of spectral a models [3], to which the Gutzwilier trace 
formula also reduces [4]. This universality has in turn provided a solid and secure ground on which 
system-specific information can be decoded by measuring deviafions of specfral correlafion func¬ 
tions from fheir universal forms, or fransifion befween fwo universalify classes. Prime examples of 
fhe former are fhe weak localization correction in Anderson Hamilfonians [5] and fhe nonuniversal 
effecf of shorf periodic orbifs (small primes) in chaofic sysfems (in fhe Riemann zeroes [6]). 

Sfudy on fhe laffer “universalify crossover", inifiafed by Dyson [7], has also come fo encom¬ 
pass a variefy of sellings, an example being fhe GUE-GOE fransifion lhaf appears in a disordered 
ring [8] and chaofic systems [9] bolh under magnetic fields. Recenf years saw applicalions of fhe 
universalify crossover in lattice QCD, in an effort lo explore fhe effecls of fhe isospin chemical 
polenlial [10] and of fhe finile lattice spacing in fhe Wilson Dirac operalor [11]. These sludies have 
revealed fhe power of fhe specfral approach in delermining fhe pion decay conslanl and fhe Wilso¬ 
nian chPT conslanfs from relalively small lattices. The aim of Ihis work is fo apply Ihis approach 
lo fhe delerminalion of low-energy conslanfs in anolher setting, namely fhe Iwo-color QCD sub¬ 
jected under fhe imaginary chemical polenlial [12] or coupled lo QED. Our novelly is lo employ 
fhe individual dislribulions of small Dirac eigenvalues [13] instead of n-level correlation functions, 
in fitting fhe lattice dafa. Practical advanfages of our mefhod will be manifested subsequenfly. 

2. chGSE-chGUE crossover 

Eel A and Bbe N /2 x N' /2 quafernion malrices, represenled by complex N xN' malrices as 
3 3 

U=h---,N/2,k=l,...,N'/2). (2.1) 

g=o^ ' g=o ^ ^ 

Here a sef of four 2x2 malrices = (I2, —io') spans fhe basis of fhe quafernion field H. Eel fhe 
malrix elemenfs belong lo A^.^^ G M and G C, so lhaf fhe malrix A is qualernion-real and B 
is nol (i.e. a generic N x N' complex malrix). We consider A^.^\ ReB^j^\ and lo be inde- 

pendenl random variables dislribuled according fo fhe Gaussian dislribulions and , 

respectively, and inlroduce an ensemble of {N + N') x {N+N') Hermilian malrices H of fhe form 

^ f n ^ V c = e-"A + (2.2) 

Y C Mn'xN' J 

Here a real parameler T plays fhe role of ficlilious time for fhe Brownian molion of fhe eigenvalues 
[7]. This ensemble enjoys fhe chiral symmefry {H, 75} = 0 wilh 75 = diag(l^, —Iat'), implying lhaf 
fhe speclrum of H consisls of A ± pairs of nonzero eigenvalues and V = |A' — A| zero eigenvalues. 
The presence of B violates fhe quafernion-realily of C and fhe selfdualily of H, tiffing fhe Ki'amers 
degeneracy of nonzero eigenvalues of H. Accordingly Ihis ensemble inlerpolales fhe fwo limiting 
cases, chiral GSE al T = 0 and chiral GUE af T —)• 00, depending on a single parameler T. 
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We consider the case in which the Kramers degeneracy is weakly broken by T <C 1. Then the 
spectral density of H in the largeW limit is identical to that of the chGSE (t = 0), i.e. Wigner’s 
semi-circle p (A) = \/4N — X'^/n. We magnify the vicinity of the origin of the A axis by introducing 
unfolded variables r,- = A,/4 with 4 = l/p(0) = 7r/\/4iV. In order to realize a nontrivial crossover 
behavior, we take the triple-scaling limit N,N' —)• oo, A, —)■ 0, T —)■ 0 while keeping the combinations 
p = ^/rjA, V =N’ —N{> 0), andr,- fixed finite. Then the j.p.d. ofN positive unfolded eigenvalues 
Pn{xi,... ,xn) is expressed as a Pfaffian of the dynamical Bessel kernel K{x,y) [14], 


S{x,y) I{x,y) 
Dix,y) S{y,x) 


Z = 


0 1 
-1 0 


1, (2.3) 


PN{xi,...,XN)=Pf(z[K{xi,Xj)]fj^^'j , K{x,y) = 

. f Jv{7i:x)yJy-i{7iy)-xJy-i{ nx)Jy{ny) _ Jy jy Vy(uy) j 


5(r,y) = 71^ 


[■ 


x^ — y2 


D{x,y) = [ dvv [ {Jy{vux)Jy{vy) —Jy{vx)Jy{vuy)} , 

2 Jo Jo 


^{x,y) = 


y/xy 


dvv^eT^^ {Jv(vx)yJv-i(vy) -xJy-i(vx)Jy(vy)}. 


Due to the recursion relation dx^Pf (^Z lK(xi,xj)]^ = {N — k+l)Pf(^[K{xi,Xjy\lj^^ , cor¬ 
relation functions of n eigenvalues are given by 

N\ f n \ 

Rn{xi,...,Xn) = dXn+l...dXNPN{xi,...,XN) =Pf [Z[K{xi,Xj)]lj^.^j . (2.4) 


3. Individual eigenvalue distributions 


The Pfaffian forms in (2.3)~(2.4) originafe from quafernion deferminanfs (Tdef) composed of 
a quafernionic kernel, [Ji^{xi,Xj)]- j, whose C-number represenfafive is fhe anfisymmefric mafrix 
Z[K{xi,Xj)]. j. Accordingly, fhe probabilify Ek{s) for an interval [0,^] fo confain exacfly k eigen¬ 
values is also given in terms of fhe Fredholm Tdef of a quafernionic infegral operafor i.e. fhe 
square roof of fhe corresponding Fredholm deferminanf of Ky (i.e. Fredholm Pfaffian of ZK^), 




(3.1) 


Here Ky denofes an infegral operafor wifh fhe dynamical Bessel kernel K{x,y) (2.3) acfing on fhe 
space of fwo-componenf L^-funcfions over fhe interval [0,^]. Firsf few EtisYs are expressed as 

^2 

/ '7 \ 

(3.2) 


Fo(^) = Del(l-^,)i/2, Ei{s)=Eo{s)^, E2{s) = Eo{s) ^ , 


3! V 8 


Fa (s) =Eo{s)U^-lnT 2 + Ta ) , E^is) =Eo{s)U^- + ^7^ + IT.n - 3 tA , 


4! V 16 4 


where T„{s) = Tr(Fi(/— denote functional fraces of fhe resolvenfs of Ky. Probabilify 

disfribufion pkA of fhe A:* smallesf positive eigenvalue is fhen given as pk{s) = 

An efficienf way of numerically evaluafing fhe Fredholm deferminanf of a frace-class operafor Ky 
acfing on L^-funcfions over an interval [0,r] is fhe Nysfrom-fype discretization [15] 


Del(l —Ks) ~def(/ —Ky), 


K, = [K{xi,Xj)y/WiWj]’”_.^^ . 


(3.3) 
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Here we employ a quadrature rule consisting of a set of points {xi} taken from the interval [0,^] and 
associated weights {w,} such that jQf{x)dx ~ YULifi^dwi- Similarly, the resolvents in (3.2) are 
approximated as Tn{s) ~ tr (Ks(/ — . For a practical purpose we choose the Gauss quadra¬ 

ture rule, i.e. sampling {x,} from the nodes of Legendre polynomials normalized to [0,s]. Previ¬ 
ously we applied the Nystrdm-type method to the dynamical Bessel kernels interpolating chGSE- 
chGUE (2.3) and chGOE-chGUE and evaluated the smallest eigenvalue distributions pi{s) [12]. 
In this work we extend our computation to the first four eigenvalues, aiming to reduce the fitting 
errors in determining the low-energy constants. We set the approximation order m to be at least 
100, and confirmed fhe sfabilify of fhe resulfs for increasing m (up fo 200 ~ 400). 

The disfribufions ,P 4 {s) for fhe v = 0 case, compufed from (2.3)~(3.3) for p < 0.70 

are exhibifed in Eig. IE. A practical advanfage of using individual eigenvalue disfribufions over 
fhe specfral densify R\ {x) = Pk{^) = S{x,x) (Eig. IR) for tiffing fhe lattice dafa is clear from 

fhe figure: fhe oscillafion of fhe laffer immediafely becomes sfrucfureless and insensifive fo fhe 
inferpolafion paramefer p due fo fhe overlapping of multiple peaks of fhe former, whereas fhe quasi- 
Gaussian shape of each peak is clearly disfinguishable and is exfremely sensitive fo p. Anofher 
advanfage specific fo fhe currenf case originafes from fhe facf fhaf P 2 k-i{^) and P 2 k{^) move in 
opposife direcfions as p is increased fo break fhe Kramers degeneracy. By combining fhe fwo besf- 
fiffing values of p for fhese fwo disfribufions, any error presenf in fhe mean level spacing A of fhe 
Dirac specfrum, which would resulf in shifting fhe unfolded dafa of 2/: — 1* and Ik'-^ eigenvalues fo 
fhe same direction, is expecfed fo be cancelled. We have confirmed fhis by generafing 10^ samples 
of crossover random mafrix ensembles wifh N = N' = 6^ and various p < 0.50 and by tiffing 
histograms of firsf four eigenvalues fo fhe analyfic resulfs. Combined values of p from fhese four 
fiffings have reproduced fhe frue inpuf values wifhin a few per mil of sysfemafic error (max. 0.5%), 
an order of magnifude closer to fhe inpuf values fhan using any single individual disfribufion. Such 
an accuracy could neifher be hoped for had we used fhe specfral densify R\ (x) for tiffing. 

4. Effective theory and low-energy constants 

The Dirac operator ^ of a QCD-like fheory wifh quarks in a pseudoreal (real) represenfafion. 



Figure 1: First four eigenvalue distributions pi,^ 4 {s) (left) for 0.04 <P<0 .70 (step 0.01, purple to red) and 
the spectral density Ri (.r) for 0.01 < p < 1.00 (step 0.01) for the chGSE (black) to chGUE (grey) crossover. 


4 











Individual eigenvalue distributions for chGSE-chGUE crossover 


Shinsuke M. Nishigaki 


such as the fundamental of Sp{2N) {SO{N)), possesses an antiunitary symmetry unlike QCD with 
quarks in a complex representation [16]: Ip commutes with ^Z* C^*), with ‘tf being the charge 
conjugation and * the complex conjugation. As = +1 (('^*)^ = — ^),D can be brought to 

a real symmetric (quaternion selfdual) matrix by a similarity transformation. Due to this property, 
the distinction between left-handed quarks and conjugated right-handed quarks is lost, leading to 
the Pauli-Giirsey extension of the flavor symmetry from S\]{Nf)l x S\]{Nf)r to S\J{2Nf) =■ G 
and its vector subgroup from S\J{Nf)v to Sp{2Nf) or SO{2Nf) =■ H. Accordingly its low-energy 
effective theory becomes a nonlinear a model on an exotic Nambu-Goldstone manifold G/H. 

Since the Dirac operator charged under the U(l) gauge field is complex, coupling QCD-like 
theories with electromagnetism or even subjecting them to the constant U(l) background breaks 
the antiunitary symmetry of ip and the Pauli-Giirsey extended flavor symmetry. In the latter case 
that is equivalent to putting on a weak imaginary chemical potential /r = //T/, its effect on the low- 
energy Lagrangian is systematically incorporated by the flavor covariantization of the derivatives 
[17]. Furthermore, if the theory is in a finite volume V =LP and the Thouless energy Ec ~ F'^/'LL?- 
is much larger than m, the path integral is dominated by the zero-mode integration (the e regime), 

Z= f dU &xp{\vLmRtivMU-VpL}F^tv{BU^BU+ M)] . (4.1) 

7su{2Ay) \2 j 

Here F is an S\J{2Nf) matrix-valued Nambu-Goldstone field, B = 03 ( 8 ) 1 ^^ ,M = ia 2 'SilNf (c^i ® ) 

for quarks in a pseudoreal (real) representation. L = (yy) /^f denotes the chiral condensate and 
F the pseudo-scalar decay constant, both measured in the chiral and zero-chemical potential limit. 
Note that the above OD a model for the case of fermions in a real representation can as well be 
derived from the random matrix ensemble ( 2 . 2 ) through the standard procedure: (i) introduce Nf 
species of complex Grassmannian {N -|-A')-vectors i//) , Yf and consider a replicated spectral de¬ 
terminant (det(A —7/)^^) = (^f where (•••) denotes averaging over A and 

B, (ii) perform Gaussian integrations over A and B, (iii) introduce a 2Nf x 2Af-matrix valued 
Hubbard-Stratonovich variable Q and open up the 4-fermi term, (iv) perform Gaussian integrations 
over xj/ and t/f, (v) take the aforementioned triple-scaling limit and denote the angular part of Q (not 
fixed by the large-A saddle point equation) as U. Then the coefficients of the mass and chemical- 
potential terms are identified as VLm = inx and 2VF^ixj = By substituting m —)• /Aoirac 

which turns the QCD partition function into the Dirac spectral determinant, the former equality 
provides the definition of unfolded Dirac eigenvalues x = Aoirac/^ due to the Banks-Casher rela¬ 
tion L = n/AV. The latter equality is used to determine F^ from the slope of the pi-p plot. 

5. Fitting Dirac spectra of SU(2)xU(l) gauge theory 

As the aim of this work is to demonstrate the validity and advantage of the method and not to 
approach the continuum, chiral, or thermodynamic limit, we chose the simplest possible setting on 
the lattice side: (i) generate 10"* samples of quenched SU(2)=Sp(2) gauge fields U^{x) on an (in¬ 
tentionally) small lattice V = 6 "*, with aplaquette action at j 3 su( 2 ) = 6 /.?su( 2 ) = 0 ~ 1-75 (step .25), 
using the standard heat-bath/overrelaxation algorithm, (ii-a) multiply the SU(2) fields on temporal 
links Uo{x) by a constant phase e'^' with pj = 0.00524 ~ .05240 (step .00524), or (ii-b) gener¬ 
ate quenched noncompact U (1) gauge fields (x) under the Coulomb gauge-fixing condition [ 1 8] 
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and multiply the SU(2) fields {x) by exp(mu(i)A^ (x)), with eu(i) = 0.0004 ~ .0024 (step .0004), 
(iii) substitute the gauge fields into an unimproved staggered Dirac operator and diagonalize. 

Due to the absence of the matrix, the antiunitary symmetries of staggered Dirac opera¬ 
tors are swapped between real and pseudoreal representations [19]. Accordingly, our case with 
SU(2)xU(l) fundamental fermions indeed corresponds to the chGSE-chGUE crossover (2.2). The 
low-energy constants are determined by the following steps: (I) fit the histogram of each of the two 
smallest Dirac eigenvalues (i.e. four counting the Kramers degeneracy) of the pure SU(2) case to 
the rescaled chGSE (p = 0) prediction pk{Xk/A)/A by varying A, (II) combine two optimal values 
of A and their variances to determine A and thus L = 7i(AV, (III) fit the histogram of each of the 
four smallest unfolded Dirac eigenvalues Xk = A/t/A of (a) SU(2)-|-p/ or (b) SU(2)xU(l) case to 
the chGSE-chGUE prediction pk{xk) by varying p, (IV) combine four optimal values of p and their 
variances to determine p and thus j2)p^ jV. 

We first observe that the four values of A obtained in the step (I) are mutually consistent, giving 
rise to combined relative errors in £ that are extremely small, ~ 0.1 % (Table 1, top). One-parameter 
fittings in the steps (I), (Ill-a), or (Ill-b) are quite satisfactory, with ;^^/dof = 0.5 ~ 1.5 for all range 
of parameters in concern (exemplified in Eig. 2, above). We also confirmed our expectation that 
the best-fitting values of p for k = 1,3 and those for ^ = 2,4 have a tendency to counter-move, 
in favor of cancelling the unfolding ambiguity due to a tiny error within A . Relative errors in p 
are considerably reduced by the combined use of four individual eigenvalue distributions (Eig. 2 
below), and are no larger than ±.018(stat)±.005(sys). Einear response of p on p/ or eu(i) 
confirmed for the SU(2)-|-p/ case (Eig. 3, left), and the pseudo-scalar decay constant at various 
values of j 8 su( 2 ) is obtained from the slopes (Table 1, middle). Eor the SU(2)xU(l) case, the 
coefficients (equivalent of F^pj) of the trBU^BU term in (4.1) divided by extrapolated to 
^u(i) 0 are summarized in Table 1, bottom. Complete lattice results, and details of analytic and 

numerical computations presented in §2 and §3 will be reported in a subsequent publication. 


Psu{l) 

0 

0.25 

0.50 

0.75 

1.00 

1.25 

1.50 

1.75 


1.310(2) 

1.255(2) 

1.199(1) 

1.139(1) 

1.070(1) 

.987(1) 

.883(1) 

.743(1) 

F^a^ 

.284(2) 

.268(2) 

.247(2) 

.226(2) 

.205(1) 

.178(1) 

.153(1) 

.115(1) 

A’^At/aV4(i) 

220 ( 2 ) 

198(2) 

186(2) 

163(1) 

145(1) 

123(1) 

99.5(8) 

68 . 0 ( 6 ) 


Table 1: Chiral condensate E from quenched SU(2) [top], pseudo-scalar decay constant F^ from SU(2)H-p/ 
[middle], and an equivalent of F^pj (divided by from SU(2)xU(l) [bottom], all in the lattice unit. 
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Figure 2: Histograms of A:* unfolded Dirac eigenvalues for Psv{ 2 ) = 0-25 and eu(i) = 0.0008, .0012, .0016 
[above, left to right] and best-fitting pk{s) {k— 1 4) from the chGSE-chGUE crossover. The p parameter 

determined for each k, their combined values [dots] and statistical errors [band] are shown below each graph. 




Figure 3: Ratios p/p/ for SU(2)H-p/ at Psvii) = 0-5 [left] and p/eu(i) for SU(2)xU(l) at Psuii) = 0-25 
[right]. 

[7] E. J. Dyson, J. Math. Phys. 3, 1191 (1962). 

[8] N. Dupuis and G. Montambaux, Phys. Rev. B 43, 14390 (1991). 

[9] K. Saito, T. Nagao, S. Muller, and P. Braun, J. Phys. A: Math. Theor. 42, 495101 (2009). 

[10] P. H. Damgaard, U. M. Heller, K. Splittorff, and B. Svetitsky, Phys. Rev. D 72, 091501 (2005). 

[11] P. H. Damgaard, K. Splittorff, and J. J. M. Verbaarschot, Phys. Rev. Lett. 105, 162002 (2010). 

[12] S. M. Nishigaki, Prog. Theor. Phys. 128, 1283 (2012); Phys. Rev. D 86, 114505 (2012). 

[13] P. H. Damgaard and S. M. Nishigaki, Phys. Rev. D 63, 045012 (2001). 

[14] P. J. Eorrester, T. Nagao, and G. Honner, Nucl. Phys. B 553, 601 (1999). 

[15] E. Bornemann, Math. Comp. 79, 871 (2010); Markov Processes Relat. Eields 16, 803 (2010). 

[16] J. J. M. Verbaarschot, Phys. Rev. Lett. 72, 2531 (1994). 

[17] J. B. Kogut et ah, Nucl. Phys. B 582, 477 (2000). 

[18] T. Blum, T. Doi, M. Hayakawa, T. Izubuchi, and N. Yamada, Phys. Rev. D 76, 114508 (2007). 

[19] M. A. Halasz and J. J. M. Verbaarschot, Phys. Rev. Lett. 74, 3920 (1995). 


7 
















































